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ABSTRACT 

New high resolution CO (3-2) interferometric map of the IR-bright interacting 
galaxy system VV 114 observed with the Submillimeter Array (SMA) reveal a 
substantial amount of warm and dense gas in the IR-bright but optically obscured 
galaxy, VV 114E, and the overlap region connecting the two nuclei. A 1.8 x 1.4 kpc 
concentration of CO (3-2) emitting gas with a total mass of 4 x IO^Mq coincides 
with the peaks of NIR, MIR, and radio continuum emission found previously by 
others, identifying the dense fuel for the AGN and/or the starburst activity there. 
Extensive CO (2-1) emission is also detected, revealing detailed distribution and 
kinematics that are consistent with the earlier CO (1-0) results. The widely 
distributed molecular gas traced in CO (2-1) and the distributed discrete peaks 
of CO (3-2) emission suggest that a spatially extended intense starbursts may 
contribute significantly to its large IR luminosity. These new observations further 
support the notion that VV 114 is approaching its final stage of merger, when 
violent central inflow of gas triggers intense starburst activity possibly boosting 
the IR luminosity above the ultraluminous threshold. 

Subject headings: galaxies: interactions — galaxies: kinematics and dynamics — 
galaxies: individual (VV 114) 



1. Introduction 

Observations of luminous and ultraluminous infrared galaxies (LIRGs/ULIRGs) are the 
key to understanding the formation and evolution of galaxies and their associated star form- 
ing environment (Sanders & Mirabel 1996; Genzel & Cesarsky 2000). Tidal interaction 
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between two or more gas-rich progenitor galaxies is largely responsible for radially trans- 
porting the gas to the central kpc, condensing the gas and triggering subsequent starburst 
activity there (Mihos & Hernquist 1996). It is now widely accepted that the elevated level of 
infrared luminosity originates from the reprocessed emission from the dust particles surround- 
ing the starburst or AGN. Based on a number of similarities, the infrared bright galaxies 
are thought to be the local analogs of the high redshift sub-mm sources discovered using 
SCUBA on JCMT during the last decade (Blain et al. 2002), and understanding the nearby 
LIRGs/ULIRGs population is an important step toward better understanding the sub-mm 
galaxy phenomenon. 

Investigation of LIRGs/ULIRGs using molecular gas at the low J-transitions (CO (1-0) 
at 2.6 mm and CO (2-1) at 1.3 mm) were carried out extensively in the past (Bryant & 
Scoville 1999; Downes & Solomon 1998). While these CO transitions were beheved to be 
a good tracer of the optically thick dense gas identifying the extent of the starburst region 
and the distribution of the fuel for such activity, more recent studies suggest that the diffuse 
inter-clump medium may dominate the CO luminosity in these low J-transition CO emission 
(Downes & Solomon 1998). Since these two J-transitions require only a small difference in 
excitation conditions, the large scale distribution of the molecular gas inferred is usually 
quite similar (Downes & Solomon 1998). The J = 3 ^ 2 transition of CO has a higher 
excitation temperature (33 K) and critical density (~ 10^ cm~^), making it a better tracer 
of the warmer and denser molecular gas of the starburst regions. In this letter, we present 
high resolution CO (3-2) and CO (2-1) interferometric maps of the IR-bright galaxy system 
VV 114 observed with the SubmiUimeter Array (SMA)^ (Ho, Moran & Lo 2004) to trace the 
distribution of the warmer dense gas and its relation to the colder, more diffuse molecular 
gas. 

W 114 is a gas rich (Mh^ = 5.1 x IO^^Mq; Yun, Scoville & Knop 1994, YSK94 
hereafter) nearby (D = 77 Mpc) interacting system with high infrared luminosity (Lir = 
4.0 X lO^^Lo; Soifer et al. 1987). The projected nuclear separation between the two optical 
galaxies (VV 114E and VV 114W) is ~ 6 kpc. Prayer et al. (1999) found a large amount of 
dust (Mdust = 1.2 X lO^M©) distributed across the galaxy with dust temperature of 20 — 25 
K. About half of the warmer dust traced in the MIR is associated with VV 114E, where both 
compact (nuclear region) and extended emission are found (Le Floc'h et al. 2002). The MIR 
spectrum also shows a sign of an AGN in VV 114E (Le Floc'h et al. 2002). Alonso-Herrero, 
Rieke & Rieke (2002) detected abundant H II regions in VV 114E and in the overlap region 
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using the narrow-band Paa images. Scoville et al. (2000) imaged the near infrared emission 
using NICMOS on-board HST and found that the highly optically obscured VV 114E (Knop 
et al. 1994) is the brighter of the two in the near infrared. Far-UV imaging using STIS 
found several hundred young star clusters in VV 114W, while no UV emission was found 
in VV 114E (Goldader et al. 2002) which suggests that most of the activity in VV 114E is 
obscured by dust and not visible in the UV emission. 

2. Observations & Results 

VV 114 was observed under good weather conditions (r23o ~ 0.08) in CO (3-2) and 
CO (2-1) using six 6-mcter diameter elements of the SMA on September 30, 2003 and on 
September 12, 2003, respectively. The receivers were tuned to 339.0 GHz (LSB) and 226.0 
GHz (LSB) which resulted in a primary beam size of 35" and 52" (FWHM) for CO (3-2) 
and CO (2-1). The array was in the "compact" configuration which includes the shortest 
and longest projected baseline of 10 and 65 meters, respectively. The largest detectable 
structure are thus 18" and 27" in the CO (3-2) and CO (2-1) lines, and the missing short 
spacing information may have substantially impacted the continuum and CO (3-2) maps 
(see below). These observations were made during the testing and commissioning phase 
of the construction of the array, and extra calibration steps were taken to scrutinize the 
array performance. Two nearby quasars, 0132—169 and 0050—094, were observed every 
25 minutes to track the instrumental gain, and absolute flux calibration was performed by 
observing Uranus. The telescope pointing drifted significantly (10-20") during the track, 
and the resulting amplitude error primarily limits the overall dynamic range of the images 
to < 20. Initial data calibration was carried out using the Caltech millimeter array software 
package MIR which is modified for SMA. The cahbrated data were imaged using NRAO 
software package AIPS (van Moorsel, Kemball & Greisen 1996). The spectrometer was 
configured with a total bandwidth of 960 MHz with a 0.8 MHz channel spacing. The final 
spectral channel maps were made by smoothing the data to a velocity resolution of about 
22 km/s. Finally, the astrometry was examined by calibrating 0050—094 (secondary) with 
0132—169 (primary), and we estimate a positional uncertainty of ~ 0.3". A summary of 
the observations is given in Table 1. Further technical descriptions of the SMA and its 
cahbration schemes are found in Ho, Moran & Lo (2004). 

Continuum images arc constructed by averaging the entire 960 MHz bandwidth of the 
image sideband data, and neither the 1.3 mm nor the 850/i.m continuum is detected with 
3(7 upper limits of 15 mJy and 150 mJy, respectively (see Table 1). The JCMT/SCUBA 
submm continuum images by Frayer et al. (1999) suggest a source extent of ~ 30" (11 kpc) in 
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Table 1. Observed and Derived Properties of VV 114 



Parameter CO(3-2) C0(2-l) CO(l-O)^ 



e (FWHM)^ 

arcsecond 2.9 x 2.2 5.0 x 3.0 6.5 x 3.7 

kpc 1.1 X 0.8 1.9 x 1.2 2.4 x 2.2 
Cont. RMS (mJy) 

Observed 50 5 - 

Theoretical 20 5 
Peak (J2000) 

R.A. 0l'^07'^47.5'^ 0l''07™47.2^ 01^07'^47.3^ 

Decl. -17°30'25.0" -17°30'25.0" -17°30'25.7" 

SAv (Jy km/sY 2275 1620 674 

Vsys (km/s) 6020 6027 6040 

Avfwzi (km/s) 330 520 520 



'^from Yun, Scoville & Knop (1994) 
'^Natural weighting 

^The integrated flux is subject to a 20% (CO (2-1)) uncertainty 
and factor of ~ 2 overestimate (CO (3-2)), mainly because of the 
uncertainties in flux and gain calibration. 
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diameter. Prom their dust model fit to the 450/xm and 850//m measurements, we estimate a 
total 1.3 mm flux density of ~ 100 mJy. Smoothing the SMA data to the angular resolution 
of the JCMT does not reveal any emission feature, and a signiflcant fraction of the dust 
emission may arise from a smooth structure larger than our shortest baseline {6 > 27"). 
Assuming the extent of the dust emission is similar to that of the CO (1-0) emission, we 
derive a 3a upper limit in 1 mm continuum of 120 mJy. 

2.1. CO (3-2) Emission 

Our new SMA CO (3-2) image, shown in Figure la, reveals a 5.0 x 2.0 kpc bar-like 
morphology with two peaks (CI and C3) occupying the two ends and a third peak C2 
located near the center of the bar. A similar structure is seen in both lower J transitions (sec 
Fig. 1 and compare the higher contours of CO (2-1) with the lower contours of CO (3-2)). 
More than 50% of the total CO (3-2) flux is detected in the compact molecular concentration 
labeled CI in the IR-bright dusty eastern galaxy (VV114E), coincident with the NIR (Doyon 
et al. 1995; Scoville et al. 2000), MIR (Le Floc'h et al. 2002) and radio continuum emission 
(Condon et al. 1990). Assuming optically thick emission and the standard galactic CO — H2 
conversion (Sanders, Scoville & Soifer 1991), its inferred molecular gas mass is 4 x 10^ M©, 
and this accounts for nearly 40% of its virial mass (10^° M© for a deconvoled source size of 
1.4 X 1.8 kpc and Af = 100 km/s), similar to what is commonly seen in the nuclear regions 
of other infrared luminous galaxies (Scoville, Yun, & Bryant 1997; Downes & Solomon 
1998). The peak C2 roughly coincides with the brightest peaks in CO (2-1) and CO (1-0), 
and the peak C3 extends to the part of the region of VV 114W where copious clusters of 
young stars have been observed previously (Alonso-Herrero, Rieke & Rieke 2002; Goldader 
et al. 2002). Despite the strong CO (3-2) emission seen in W 114E, the main disk of 
the optically bright western galaxy (VV 114W) shows little emission, possibly resolved out 
by the interferometer. The nuclear region of VV 1 14E exhibits higher line intensity ratios 
(r2i = l2i/ho and r^i = /32//10) than the rest of the system, suggesting the presence of 
signiflcantly warmer, denser gas. The inference of high gas temperature is further supported 
by the observed peak brightness temperature of AT ~ 10 K averaged over the 1 kpc beam 
area (i.e., T ~ 100 K for a filling factor / ~ 0.1). 

The absence of baselines shorter than 10 meters and the resulting limited sensitivity 
to structures > 18" in extent may account for much of the extended structure missing in 
the CO (3-2) image. For spectral line observations, this spatial filtering is not as severe 
as one naively calculates since coherent bulk motions such as galactic rotation substantially 
reduces the effective angular sizes of the emitting regions in individual velocity channels. 
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Multi-isotopic, multi-transition analysis of molecular tracers have indicated a multi-phase 
medium with a significant contribution by sub-thermally excited diffuse gas in the turbulent 
nuclear starburst regions and thus a lower conversion factor between the CO luminosity and 
the total molecular gas mass (see Aalto ct al. 1995; Downcs & Solomon 1998). Therefore, 
some of the observed morphological differences from the lower J transitions (which are more 
easily excited) may reflect real differences in the physical properties of gas. A more definitive 
characterization of the distribution and physical properties of the warm, dense gas traced in 
CO (3-2) will require obtaining a single dish measurement in the future. The total observed 
line flux of 2275 Jy km/s translates to a total molecular gas mass of 2 x 10^° Mq, which is 
about 40% of the total gas mass inferred from the CO (1-0) measurement by (YSK94), and 
this demonstrates clearly both the large gas mass and the highly concentrated nature of the 
dense (n > 10^ cm~^) molecular gas in VV 114. 

The clumpy spatial distribution of the high density gas makes the CO (3-2) kinematics 
appear complex (see Figure 2). A rotation-like velocity gradient is seen from the southwest- 
northeast with a noticeably steeper gradient near the C2. The direction of the overall velocity 
gradient is consistent with that in CO (1-0) (YSK94) and in CO (2-1) (Fig. Id). A large 
amount of extended emission may be missing by our observations, and the derived velocity 
field should be interpreted with a caution. Both the excitation requirement for the CO (3-2) 
emission and the spatial filtering by the interferometer may be isolating the kinematics of the 
densest gas, possibly in the deepest parts of the evolving potential. The total hnewidth at 
zero intensity is 330 km/s, much narrower than those traced in CO (1-0) and CO (2-1), and 
the linewidth in VV 114E alone spans over 100 km/s. The peak CI is most hkely associated 
with W 114E, and we thus infer the systemic velocity of 6052 km/s for VV 114E. 



2.2. CO(2-l) Emission 

The velocity integrated CO (2-1) map is shown in Figure lb. The CO (2-1) emission 
extends widely across the two disks comprising W 114, showing a high degree of similarity 
to the CO (1-0) (YSK94; also Fig. Ic), including both the bar-like structure 5x3 kpc in 
extent dominating the central region and a long (3 kpc) molecular tail extending out from the 
southeastern edge with no obvious optical counterpart. The higher angular resolution of the 
CO (2-1) map makes subtle features appear more enhanced, but the overall distribution is 
very similar to the CO (1-0) map. The excellent correspondence of the two is not surprising 
since both transitions are easily excited by the physical conditions typical of cold, dense 
(T = 10 - 20 K, n ~ 10^ cm'^) gas. The global hue intensity ratio r2i = hi/ ho = 2.4 ± 0.7 
is much smaller than the expected value of 4 for fully thermalized, optically thick emission 
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from a single phase gas. This is a further indication that sub-thermal excitation is an 
important concern for CO transitions in this object. 

The gas kinematics traced in C0(2-l) and CO(l-O) lines are also nearly identical. The 
total hnewidth of the CO (2-1) emission is 520 km/s which is comparable to the linewidth of 
the CO (1-0) emission. By taking the median velocity of the emission channels, the CO (2-1) 
systemic velocity of VV 114 as a whole is estimated at 6027 km/s. A steep velocity gradient 
is seen in both CO (2-1) and CO (1-0) where the southeastern molecular tail connects to 
the main component of molecular gas in VV 114, but a smaller gradient is seen along the 3 
kpc length farther south. This may suggest that the molecular tail is tidally driven, and its 
kinematics are dominated by the angular momentum of the part of the galaxy from which 
it was initially perturbed (see below). The large velocity dispersion is seen near C2 in both 
CO (2-1) and CO (1-0) emission, possibly marking the new dynamical center of the merger 
(YSK94). Alternatively, the large velocity dispersion may imply a tidally induced radial 
motion seen along the line of sight. 



3. Discussion 

While several hundred young star clusters have been identified in W 114W (Goldader 
et al. 2002), most of the current vigorous activity is seen in the highly obscured galaxy, 
VV 114E, as identified by its infrared, radio and sub-mm continuum emission. The origin 
of the intense IR emission, starburst or AGN or both, in VV 114 is still subject to debate. 
While Le Floc'h et al. (2002) suggests that 40% of the MIR emission is associated with the 
compact nuclear AGN in VV 114E, different diagnostics (i.e. NIR spectroscopy, Doyon 
et al. 1995) suggest that the IR emission is mostly dominated by young stars. The three 
massive, compact gas complexes identified by our CO (3-2) observations mark the locations 
of the densest concentrations of gas and the sites of current activities fueled by this gas. The 
dominant peak CI in W 114E appears to mark the location of the most intense current 
activity. The distributed nature of CO (3-2) emission and the presence of several massive gas 
concentrations provide the clearest evidence yet that distributed intense starburst activity 
may provide a significant fraction of the large IR luminosity in W 114. 

Our new CO (3-2) and CO (2-1) observations strengthen the earlier suggestion that 
VV 114 is a late stage merger ready to undergo more massive bursts of star formation in 
the near future. The widely extended CO (2-1) and CO (1-0) emission across the two 
galaxies implies that the cold and dense gas have begun to decouple from the gravitational 
potential of the host galaxies, and such decoupling only occurs in the well advanced stages 
(t ~ 5 X 10^ years since the initial collision) before the coalescence results in a massive 
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elliptical galaxy (e.g., Mihos & Hernquist 1996). The three discrete peaks traced in CO (3- 
2) emission identify the massive and compact molecular gas concentrations condensing out 
of the inflowing gas. By comparing the observed distribution of gas and optical morphology 
to those predicted in simulations (Mihos & Hernquist 1996; lono, Yun & Mihos 2004), we 
estimate that the final coalescence may occur within the next ~ 10^ years. The CO (1-0) and 
CO (2-1) gas distribution of VV 114 is very similar to what is seen in the IR-bright galaxy 
NGC 6090 (Bryant & Scoville 1999; Wang et al. 2004). There the widely extended gas is 
seen over the edge-on and face-on system separated by 3.4 kpc, where not only the CO (1-0) 
and CO (2-1) but also the CO (3-2) gas peaks in the overlap region (Wang et al. 2004). 
Thus, both VV 114 and NGC 6090 appears to be undergoing a similar transition where the 
bulk of the gas is in the process of decoupling from the host galaxies and funneling toward 
the dynamical center of the new combined potential, possibly resulting in the ultraluminous 
{LjR > 10^^ Lq) phase, following the scenario advocated by Yun, Scoville & Knop (1994). 
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a (2000) 



Fig. 1.— (a) CO (3-2), (b) CO (2-1), and (c) CO (1-0) (YSK94) maps overlaid over R-band 
optical image. The contour levels represent 4,6,8,10,12,14,16 xa (40 Jy km/s) for CO (3-2) 
(e = 2.9" X 2.2"), 6,8,12,16,20,30,40,50,60 xa (5.0 Jy km/s) for CO (2-1) {6 = 5.0" x 3.0"), 
and 4,6,8,10,15,20,25,30 xa (3.5 Jy km/s) for CO (1-0) {9 = 6.5" x 3.7"). The three resolved 
CO (3-2) molecular complexes are labeled CI, C2 and C3 from east to west, where CI is the 
most dominant of the three (~ 50% of the total flux). Each of the CO images are centered 
at the phase center of each observations. No primary beam correction is applied since the 
source size is smaller than the primary beam even at 330 GHz. (d) The CO (2-1) velocity 
distribution contours (Ml) (units in km/s) overlaid over the velocity dispersion (M2) map 
in gray scale which ranges from 10 to 110 km/s. The velocity contours are equally spaced in 
50 km/s intervals. This figure was made using the WIP software package (Morgan 1995). 
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1^07^48.0^ 47.5^ 47.0^ 46,5^ 

Fig. 2. — The CO (3-2) velocity distribution is shown in gray scale (in km/s) overlaid with 
the velocity contours of the CO (2-1) emission. The contours range from 5900 to 6100 km/s 
in steps of 50 km/s. 



